Inducing spin-dependent tunneling to probe magnetic correlations in optical lattices 
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We suggest a simple experimental method for probing antiferromagnetic spin correlations of two- 
component Fermi gases in optical lattices. The method relies on a spin selective Raman transition to 
excite atoms of one spin species to their first excited vibrational mode where the tunneling is large. 
The resulting difference in the tunneling dynamics of the two spin species can then be exploited, 
to reveal the spin correlations by measuring the number of doubly occupied lattice sites at a later 
time. We perform quantum Monte Carlo simulations of the spin system and solve the optical lattice 
dynamics numerically to show how the timed probe can be used to identify antiferromagnetic spin 
correlations in optical lattices. 

PACS numbers: 67.85.d, 73.21. Cd, 74.25.q, 75.50.Ee 
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I. INTRODUCTION 

Atoms in optical lattices provide a powerful test bed 
for probing the properties of quantum systems in peri- 
odic potentials. Several fundamental results have been 
reported including the observation of Mott physics [IJ, 
fermionic pairing [2 , and the fermionization of bosons 
in one dimension [3 . Promising new research directions 
are emerging with the creation of Dirac points in honey- 
comb lattices ^ and the achievement of single site res- 
olution for probing and manipulating two dimensional 
(2D) lattices [H [6]. A major objective in the study of 
atoms in optical lattices is to unravel the physics of quan- 
tum magnetism, and for this purpose it is crucial to be 
able to measure the spin correlations of the atoms effi- 
ciently. Several methods have been suggested including 
Bragg scattering [T, time-of- flight measurements [SfllQj, 
and light polarization experiments [11 . A complication 
for these methods is, however, that the experimental sig- 
nal is small, and that one has to average over many ex- 
perimental realizations. 

To address this, we recently suggested a super lattice 
method which yields large macroscopic signatures of the 
spin correlations and allows for the detection of both 
short- and long-range correlations [12]. It is based on 
spin-dependent transport using a spin-selective Raman 
excitation followed by the application of an additional 
moving superlattice potential. The spin-correlations can 
then be extracted by measuring the number of doubly 
occupied lattice sites. 

In this paper we propose a similar but experimentally 
simpler method, which replaces the moving superlattice 
by a precisely timed period of free tunneling dynamics of 
the atoms followed by a measurement of the double occu- 
pancy. We analyze the spin dependence of the tunneling 
dynamics and show how this can be used to reveal the 
spin correlations in the initial state. The main result of 
our paper is the demonstration that timed measurements 
of the probability for double occupancy of the lattice sites 
yields distinct signatures of the magnetic correlations. 



Importantly, the method we propose here does not re- 
quire any spin dependence of the lattice potential and 
only utilizes only the optical lattice potential which is 
already present. This is a major advantage from an ex- 
perimental perspective. 



II. MODEL AND METHOD 

We consider a gas of fermionic atoms in an optical 
lattice with equal populations in two internal states |cr) 
{a =^, ^), which we refer to as spin in the following. For 
strong repulsion and half-filling, the gas is in the Mott 
phase at low temperature T which is well described by 
the Heisenberg model [13j 



H = J 



{l,m) 



{l + A)sts^J. (1) 



Here, s^ is the spin- 1/2 operator for atoms at site /, 
and (l^m) denotes neighboring pairs. The interaction 
strength J can exhibit a small anisotropy, quantified 
by the anisotropy parameter A. A major incentive for 
studying atoms in optical lattices is to observe the phase 
transition to an antiferromagnetic (AF) ordered state 
which occurs at T = T/v — 0.945 J in 3D for the isotropic 
(A = 0) Heisenberg model [9]. To achieve this, one needs 
a sensitive experimental probe of the spin correlations. 

The key idea of the present paper is to access the spin 
correlation functions by inducing a difference in the tun- 
neling dynamics of the two different spin states 1 1) and 
I D as shown in the sketch of Figure Ilia). To this end, 
we excite all 1 1) atoms to their ffist excited vibrational 
state in the lattice sites while the | l) atoms remain in 
the ground state [14 . This can be done experimentally 
by set of spin selective Raman transitions on the vibra- 
tional sideband of an auxiliary state. Then we simply al- 
low the system to evolve for some time and subsequently 
measure the number of sites which are doubly occupied 
Nd. Since the tunneling to neighboring sites is larger 



for the first excited state tiian for tiie ground state, tiie 
I D atoms essentially remains at the initial lattice site as 
shown in Figure flla). In practice, the number of dou- 
bly occupied sites can be counted by merging two atoms 
residing on the same lattice site to a molecule by ap- 
plication of a photo-association pulse or via a Feshbach 
resonance. [T, TF 

In our theoretical modeling the lattice potential is held 
constant during the procedure and we assume that the 
interaction is tuned to a sufficiently weak value, for in- 
stance by using a Feshbach resonance [16 , so that the 
system can be considered non-interacting during the tun- 
neling dynamics. The three spatial components of the 
atomic wave function then evolve independently. Assum- 
ing the motion in the y- and z-directions to be frozen 
by strong lattice potentials, it is sufficient to consider 
the dynamics along the x-direction, where the atoms are 
trapped by a periodic potential A{x^t) = Asm^{7rx/d) 
with d the lattice constant. With this potential the single 
particle Bloch eigenenergies £^[ and corresponding eigen- 
functions \ipl.) are indexed by a vibrational state label 
/ = 0,1,2,... and a Bloch wavevector k. The initial 
states are taken to be the Wannier functions given by 
\Wj^) = -^ Zl/e ^~*^"^^l'^fc)' which describe a localized 
state of the n'th vibrational mode centered around lat- 
tice site j [17] . 

Because the Wannier functions are not eigenfunctions 
of the non-interacting lattice Hamiltonian, they evolve 
over time. To describe this we introduce the transition 
matrix U{t) = e''^' = Y^ijn ^If {-t)d]^di^, where a]^ is 
the creation operator for the Wannier state |W.^) and we 
use units for which h = 1. The time evolution operator 
tj{t) is quadratic in a and d) and diagonal in the vibra- 
tional quantum number n since we assume the system 
to be non-interacting during the evolution. The matrix 
elements of this transition matrix read 



Eq. (|3| reduces to 



u^\t) = {w-\wr{t)) = ^Y.^"'^'~ 



■l)d. 



-iEV:t 



(2) 



where \Wl^{t)) represents the initial Wannier function 
\W{^) evolved for a time t. We are mainly interested 
in the probability for a site i to be doubly occupied 
Pd{i) = (^it(^)^ii(^))- 0^^ proposed scheme is based on 
the fact that this probability can be expressed in terms 
of spin correlation functions of the initial state and the 
overlap matrix elements Ui^ [12 . Since the ||) atoms 
remain in the vibrational ground state n = while the 
It) atoms are transferred to the first excited state n = 1, 

we write [/• • = Uj, and U>- = Uf,. Then we obtain 
''J ''J f'j ''J 



pa{i) =Y.Hit)\"pfm'' {\ 



{sUf}) 



*j*i 



ji 



-Hululufpf^is^si). 

3^1 



P'^^) = j-Y.Hit) 



m 



(4) 



Equations ([3| and Q demonstrate the direct link be- 
tween the observed double occupancy and the spin cor- 
relation function (s|sf) describing AF ordering in the 
initial state. 



III. RESULTS 

We now present our numerical results. The spin corre- 
lation functions (si-Sj) for the initial state determined by 
the 3D Heisenberg model ([l]) are obtained from full quan- 
tum Monte Carlo (QMC) simulations using the stochastic 
series expansion method [18] with directed-loop updates 
[T9^ . The subsequent time evolution after the 1 1) atoms 
are excited to the first vibrational state is then calcu- 
lated by numerical diagonalization of the non-interacting 
Hamiltonian taking the motion to be along a single direc- 
tion. As an example, we will here present results for a ID 
lattice with A^ = 32 sites and a lattice height A = lOEji 
where Er = h^n'^ /2Md^ is the recoil energy with M 
the atom mass. The numerical diagonalization is per- 
formed using a lattice unit cell consisting of 1028 real- 
space points. 



A. Single atom dynamics 

We first consider the time evolution of a single atom. 
Figurewmb) andmc) depicts the probability \UQAt)\^ of 
a |<j) atom which siarts at t = at site z = 0, to be at site 
j at time t. We see that for the parameters chosen, the 
I \) atoms do not move but remain at their initial position 
(Fig. [lib)) whereas the |t) atoms move (Fig. Illc)). This 
makes this tunneling dynamics well suited for measuring 
AF ordering. Note that at certain times the wavefunc- 
tion of the 1 1) atoms mostly samples either even or odd 
lattice sites in a way similar to a random walk in ID. In 
Fig. ^b) we plot the tunneling time of the | t) atoms, 
t^, defined as the first zero of 9^|/7QQ(t)p, as a function 
of the lattice depth. According to the WKB approxima- 
tion this time depends exponentially on the lattice height, 
and we indeed find t^ ex exp(0.4A/£^i^) for deep lattices. 
Thus, by adjusting the lattice depth one can tune the 
time scale of the probe. It is, however, important that 
the lattice depth is sufficiently deep so that only the 1 1) 
atoms move. This is illustrated in Fig. |2] (right) where we 



(3) plot Pel 



^ii (^^)p5 which represents the proba- 



In the limit where the | \) atoms do not move | Ui 



■^12^ A. 



bility that a | \) has moved after the 1 1) tunneling time t^ . 
As expected, the error decreases rapidly with increasing 
lattice depth. 
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FIG. 1: (a) A cartoon depicting the tunneling dynamics for 
the two spin species, (b-c) The actual tunneling dynamics 
for the two lowest lying vibrational Wannier states of a one- 
dimensional lattice with lattice height A — IOEr. We show 
the norm square of the wavefunction projected onto Wan- 
nier states centered around each lattice site. There is a clear 
difference in the tunneling dynamics between (b) the almost 
stationary vibrational ground state (| \.) atom) and (c) the 
first excited state (| t) atom) which tunnels to its neighbor- 
ing sites. The dashed line indicates the first minimum of the 
probability for the 1 1) atom to remain at its initial position. 
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FIG. 2: (Color online) (a) the tunneling time versus lattice 
depth A of the first excited vibrational state which sets the 
time scale of the probe, (b) the error probability as a function 
of lattice depth A. 



B. Many-body dynamics and AF ordering 

We now analyze the evolution of the double occupancy, 
Pd starting from a correlated initial state determined by 
the 3D Heisenberg model. This is shown in Fig. [3] 

We focus first on the non-interacting case correspond- 
ing to T = oo which is plotted as the solid line forming 
the boundary to the grey area. For this case, the spins 
are uncorrelated and (s^s^) = ^Smw Since the vibra- 
tional ground state is essentially immobile \UJ'j\'^ ^ 6ij 
for all the relevant times, we find Pd ^ ji^ ^ irn f+\\'2\ 



\uUtW 



be larger than 1/4 and will in general be lower because 
of the non-vanishing probability that the 1 1) atoms re- 
turn to their initial states. The oscillations in Pd for 
uncorrelated atoms are therefore directly linked to the 
oscillations in \Uj^{t)\'^ which are shown in Fig. hi Since 
the tunneling time f^ gives the time for the first minimum 
of \Uj-{t)\'^^ it corresponds to the first maximum of P^. 
At this minimum, \Uj-{t)\'^ c^ and the value of Pd at 
the first maximum is therefore very close to 1/4, which 
is also the long time uncorrelated value. 

We now turn to Pd for finite temperature also plotted 
in Fig. [3J We have chosen a slightly anisotropic Heisen- 
berg system with A = 0.01, which forces the staggered 
magnetization to be along the z-axis. The AF correla- 
tions are clearly visible as increased oscillations of Pd as 
a function of time which is the main result of the present 
paper. In particular, note that the probability Pd is larger 
than 1/4 for certain times which is a clear sign of AF or- 
dering. To analyze this effect further, we consider the 
high temperature regime J/T <C 1 with /c^ = 1, where 



i^tsj) 



1. 



J 

i&r 



Si.j) + OiJ^T^) 



(5) 



with S(^ij^ the Kronecker delta function connecting neigh- 
bor sites. Using this in Eq. Q we get 

P,{i)=^-{l-\Ul{t)\^) + ^\ul^,^{t)\^ (6) 

for J/T <C 1. Equation ^ clearly illustrates how AF 
correlations can increase Pd] for example as seen from 
Fig. [l] at t = t^, we have that \Ul{t^)\'^ c^ and 

\Ul±i{t^)? - 1/4 such that Pd{i) = 1(1 + ^). As the 
temperature is lowered, the oscillations of Pd becomes 
larger refiecting the increasing AF correlations. 
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FIG. 3: (Color online) The double occupancy probability Pd 
as a function time t at different temperatures T for the slightly 
anisotropic Heisenberg system (A = 1.01) and a lattice depth 
A = IOEr. 




FIG. 4: (Color online) The difference in the double occupancy 
probability APd between the interacting and non-interacting 
cases as a function of time t and temperature T. 
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FIG. 5: (Color online) The /c-space dependence of the norm 



square of the transfer matrix elements Uq- for the first excited 
state shown at selected times. Note how the timed tunneling 
dynamics allows to probe the spin structures other than k — 
TT, which is relevant for AF ordering, like the two examples 
shown here which mainly probe k = 7r/2 and k = 7r/3 ordered 
states. 



We see that for the oscillations at later times, the dif- 
ference increases significantly for T < T/v ~ 0.945 J when 
long-range order sets in, although the absolute size of 
the difference is rather small. Note also the two dips at 
t ^ 260h/Eii and t ^ bOOh/En which shows up as long- 
range correlations set in at T = Tn- From Fig. |4] one can 
also see that short-range correlations leave a signature at 
T > T/v, for example yielding an AF signal at t ^ t^ of 
APd{T = 2J) = 0.0195 ^ 0.08P(T = oo). 



IV. CONCLUSIONS AND DISCUSSIONS 

In summary, we have demonstrated how to detect anti- 
ferromagnetic spin correlations in optical lattices by uti- 
lizing the difference in the tunneling dynamics of the 
two lowest lying vibrational states of the atomic wave- 
functions. Measurement of the number of doubly occu- 
pied sites yields a macroscopic signature of the magnetic 
correlations which is proportional to the number of lat- 
tice sites. Below the critical temperature T/v the method 
has large signatures AP^ ^ 0.1 as a result of the spin 
correlations. 

The method presented here is particularly relevant for 
measuring short range correlations as it provides an easily 
accessible experimental method, which requires a mini- 
mal change of the experimental setups already present. 
For measuring long range correlations the method is less 
suited as the spreading of the wave function for long 
times introduces an uncertainty in the correlation func- 
tion measured. Such long correlation can thus be mea- 
sured more precisely with different schemes, for instance 
relying on superlattice potentials [12]. 

Here we have mainly focussed on the detection of anti- 
ferromagnetic correlations, but the method may also be 
used to probe other spin ordered states. In Fig. [5] we 
show the /c-space representation of \UJq\'^ for three dif- 
ferent selected times. As seen, by selecting appropriate 
times a magnetic ordering with periodicities of, for ex- 
ample, k = 7r/2 and k = 7r/3, can also be probed. 



We note that the oscillations of P^ increase signifi- 
cantly when T < Tn where the system is in the broken 
symmetry phase. This is particularly visible for the later 
times which indicates the onset of long-range order. This 
effect is illustrated further in Fig. [4j where we plot the 
difference AP^ between the probabilities for finite and 
infinite temperatures. 
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